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On the Bilcclation of certain integrals occurring in the theory 
of molecules, especially three-centre and four-centre integrals 


By Stic O. Lunpevist and PErR-OLtov LéwpiIn 


By using a method involving expansion in spherical harmonics, it is shown that a general 
formula can be deduced for the potential integrals, which occur in the theory of molecules 
and crystals. Explicit formule are given for the two-centre, three-centre and four-centre 
integrals of atomic S-functions. The method can be applied to analytical wave functions, 
but the formulz are more convenient for computations by means of numerical integration. 


I. Introduction 


In the theory of molecules and crystals the evaluation of certain integrals 
over products of atomic wave functions, particularly exchange integrals, three- 
centre and four-centre integrals, presents considerable difficulties. ‘enone of 
these types have been treated by many authors, [1]—{30], but, as far as we 
know, any general way of investigating them has hitherto not been developed. 

For analytical wave functions of the hydrogen-like type, all two-centre integrals 
can conveniently be évaluated by using elliptical coordinates, [2]—[17]. Some 
three-centre integrals have also been treated in this way, [14]—{17], but the 
results were here confined to some special configurations of the nuclei. A general 
formula for the three-centre integrals has not been given', and the four-centre 
integrals seem to have been neglected in a large part of the literature.” 

The methods used here are certain expansions in spherical harmonics, first 
introduced by CoonipGE [9] and later used extensively by LANDsHoFF [22] and 
by Lown [19]. The purpose of this note is to show that these methods can 
easily be generalized to include all the potential integrals of the type considered 
here and that general formule can be deduced even for the three-centre and 
the four-centre integrals. Explicit formule will be given for S-functions.® 


1 For some special cases see [9], [18] and [19]. 

2 An approximate formula has been given by SKLAR [28]; cf. also [29] and [30]. 

3 The formula for the exchange integral given here, (24), was first deduced by LaNpsHorr. 
The formule (23)—(28)- were originally derived by the present authors independently of each 
other. Special cases of the formule for the three-centre integrals were already deduced by 
CooLipGE in treating a molecular problem. The necessity of taking the terms corresponding 
to the three- and four-centre integrals of the type considered here into account in calcula- 
tions of the cohesive energy of ionic crystals was shown by FrROMAN and LuNp@vistT Ark. f. 
fysik 2, N:o 40 (1950). 
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II. The general potential formula 


Let us consider a potential integral of the general type 
ii, Q, (1) Oe a ae (1) 
T12 


where (2, (r,) and (rz) are arbitrary density functions, dt, and dt, are the 
volume elements, and 7,. is the distance between the points 1 and 2. In order 
to evaluate (1) we will use the expansion of 1/r;, in spherical harmonics. If 
r, and r, are vectors from the same origin, « is the angle between them, and 
vy; and 7, denote the smallest and largest of the quantities r, and 72, we have 
the expansion 


= >) Fa Pi (cosa). (2) 


Introducing a fixed system of spherical polar coordinates (7, 3, y), and using 
the addition theorem for the Legendre polynomials, we get 


2h+1(h— my)! 
Spe o, = kpm P™ tam @ n= 
im(D, 9) ~ kim Pp (008.8) "9, fam = |/ EG (4) 
They are orthonormalized and satisfy the relation 
27 7 
[de | Sim Sn m sin 9d9 = Onn Om'm- (5) 
0 0 


Putting the expansion (3) into the integral (1), we obtain directly the general 
formula 


Ri 
+) 2P (4) 1 das jo8 (7) "9% ar} = 


0 i 
~~ athe we 20 . (6) 
~ 2 2 aH pif Male my drs [O58 (ra) dry + 
se (Fyre dry pa, (19) ° 7 12dr,l 
0 0 
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where the last term has been transformed by using Dirichlet’s formula for double 
integrals. The quantities Qr, , (N = 1,2) are here given by the integral 


QR (7 - fag | ox rd, p) Sam (9, y) sin dd 9, (7) 


and they can therefore be interpreted as the coefficients of the expansion of 
Qy into spherical harmonics. Formula (6) forms the basis of our method of 
calculation. 


III. Applications to two-centre, three-centre, and four-centre integrals of 
S-functions 
In the theory of molecules and crystals, a fundamental integral is of the type 
I III IV 


ibacike fee yea rel ee (8) 


Ty2 


where the indices a, b, c, d denote the different nuclei involved, and the Roman 
numerals indicate wave functions of various types. This integral is now a 
special case of the potential integral (1) with the density functions 


i IL ny, 


Q, = Ve vie Q, = Je5 ae (9) 


and (8) can therefore be evaluated according to the general formula (6). For 
the sake of simplicity we will here confine ourselves to considering only S- 
functions, but we note that the method can be extended to include also P-func- 
tions and higher functions, even if, in the latter case, the calculations will be 
more complicated. 

Let us choose the point a as the origin of our coordinate system (r, 0, ¢), 
and let the vector rac represent the direction of the polar axis. Let the plane 
g=0 be determined by the vectors race and rap, and let further ® be the 
angle between this plane and the plane determined by the vectors rac and aa, 
1.e. ga = ®. The S-function ¥Y, is spherically symmetrical around the nucleus 
y (v=a, b,c, d) and may here be defined by the expression 


W,(1) = Trek (711). (10) 


According to (7) and (9) we obtain directly 


1 ir ou 
Or, (a) = ia fan [amin a) 9b (715) Sam (01, V1) sn 7 dD, 
0 0 
(11) 
| ra “1 IV 
257 ("2) = fix | daft (720) ga (T2a) Sim (Pg, Pg) Sin Dy d Ip. 
0 6 
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In order to evaluate the integrals (11), we will now expand the functions g, 
into spherical harmonics around the point a in the system (7, 3, vy), which for 
S-functions can be carried out in a simplified way. If #wa, denotes the angle 
between the vectors ry and r,, we obtain at first the following expansion into 
Legendre polynomials 


oo 


(ry ») an (Tav | rn) [2h (cos ya»), N — Ae Z 
h=0 
(12) 
2h+1 
an (Tav | tw) = aah tale (rv») Pr (cos nay) Sin Ovary d Onay- 
0 
For » = 6, c, and d we then get the expansions into spherical harmonics 
iat Soe Cie ra 
go (710) a 2 oh +1 nh (rab |11) Sim (b) Sim (91; Pr); 
IL % TIT 
ge (2c) = >) &n (Tac| 72) Pn (cos Bg) (13) 
h—=0 


IV ny 


oo +h 
ga (724) ei pe re awe) tn ( (raa| 2) Sim (4) Sim (92, G2); 


where, in the cases y= 6b and »=d, we have used the addition theorem for 
the Legendre polynomials. Putting these expansions into the integrals (11) and 
using the relation (5), we obtain the following explicit expressions 


( Bane x I II fa 
| Dhe1e® (cos Oy ae) Ja (74) on ( (rav | 73); bx i: 
m J ky 1m : I 
20 (1) = Shot j m0 Ya (711) xn (Tav| 11), bec (ie 
kp, I II 
[ah 4 1 Ome 0x0 4a (Ps) da (12), reas 
and 
3 knm (L—)! yn peat nt : 
fel 2 h ey : i m)! 23 (cos Broad) é a] (Tada | 1) a) (Tae | To) dx i 
Q-m =) we T1V 
2,h (72) 9 j a l Om0 Kh (fac\ fs), d =¢ (15) 
kim It IV 
DF) at Omo OD (rac| 12) Ja (79). d=a 
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In addition to (12), we have here used the notations 


a 


TI, IV 2h + ae ov 


an (Tac| 12) = : Je (T2e) Ge (r2c) Pn (Cos Dy) sin Vy d Do, 


0 
(16) 
III “IIL 
am", (Tac lr) = aati Je (r2c) P™ (cos By) P™ (cos Bg) sin Fy d Vp. 
j 


The main problem is now the calculation of the radial functions « (r) defined 
by (12) and (16). Let us first treat the integral in the right-hand member of 
(12). Following Léwory, [19], p.19, we will introduce the new integration 
variable ry, = R, which gives 


R? = ri, + ry — 2 ry rar C08 Ona, 


(17) 
: R 
sm 0yard Ovary = —— aR, 
TN Tas 
and 
Tayt? 
Qh+1 r+ ra, — R 
an (Ya |r) = a ite [o (R) Pr (a=) “RAR. (18) 
For h =0, 1, ... we obtain in particular 
Tach 
1 
a0 (Tar |r) = orn | PR RGR, 
Iray—7 | 
(19) 
Fg tt Taptt 
3 APY t ‘ 
Hy (rav|1) = 4p re Gs = ri.) | gv» (R) Rd R — | gv (R) Fans, 
lav? tao 


More complete formulae may be found in Lowp1n, [19], p. 20-21. This means 
that the radial functions «(7a,|r) can be easily calculated from the functions 
g,(r),. whether the latter are given analytically or numerically. The formulae 
(19) are, in fact, particularly suitable for computations using numerical integra- 
tions. III, IV 


We note that the corresponding formulae for «(rac|7) may be obtained by 
I IV 


putting the product g- ge instead of g, into the relations (18) and (19). 
Finally, we have to treat the integral (16”). Putting the expansion (13”) 
into (16), we get at once the formula 
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TIE 2h : i Ir 
a (rac| re) = pti AM 4, Ot ( (Tac| 72) (20) 

with 
Ane = [Pew u) Py (u) du, (21) 


The coefficients A”, may be found by using the elementary rules in the theory 
of the spherical harmonics; the first few of them are 


2 
ay 52 ea Ont eg oT 4 KALE 
Hee 4 [ictal oh ei ola ai (22) 
i= Di Qiao ere 
21 (Ie A oa 1 


Other quantities of this type can easily be calculated by means of the formulae 
for products of spherical harmonics! given by Cooxrpce, [9], p. 204—206. Hence, 
all the functions «(r) can be found by means of the elementary functions given 
by (18) and (19). 

Combining (14) and (15) with the fundamental formula (6), we can now 
obtain explicit expressions for all two-centre, three-centre, and four-centre inte- 
grals of the type (8) involving only S-functions: 


Two-centre integrals: 
t 


I IL WIIV ILIV “Ty 1 
(dae €) = | v0 (relra) dre | Ga (12) Ja (72) 724% + 


0 0 
(23) 
III, IV a II 
a | ao (rae |11) dn | (72) ga (72) 72d 7g. 
0 n 
I IL I1IV 00 l { my IL r| III 
(acc a)= i Qh fs 1? | Ga (1) &n (Tac| 73) rt dr, | ae(re) on (Tac| 12) ht? dry + 
h=0 « 
0 0 
(24) 


rere) iit i, - 
ze [sor tn ( (rae| 1) 7 ie ty [oe (72) on ( (rac| 2) ridin) 
0 0 
1 For the case 1=h, sce also H. Betue, Ann. d. Phys. 3, 133 (1929), particularly p. 188—192, 


| 


: 
i 
( 
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r 


I IWIM1IV Iv Ur I iat 
(aac a) =| gulrs)ao(raclrs) rds f da (12) Gora) dre + 
0 0 (25) 
co 
‘Iv Ir 1 
a [ a(n) (rac | 11) ridin | (re) dn (re) 72d 
0 at 
Three-centre integrals: 


I ILIV eo 
(aac d) = 235 {Pe (cos Bead): 
0 rs 
IIL IV I IL 
ff catredlny zn (taalra) adr | do("s) dn (re) r2dry a 
0 6 
“TH ri e (26) 
+ fea (Tac\1) ma (raa| 11) fdr | (12) ga (T2) Te ar}: 
r 
= ie 
abc a) : 
rat ve 
co TT 
“I I TIL IV 
a Ga (71) %n (rav|71) adn [ 2) (rac| 72) Ga (Te) ret? dra + 
0 6 
1 (27) 
“TI IV IL 1 
ss [ ca(reclny ) 9a (74) 7 dn [ ool) on (Yad | 12) 1 Eee ey 
0 6 


Four-centre integrals: (€9 = 1, &m = 


= 1) 
a | Pi’ (cos Bac) P7” (cos Paac) cosm® - 


BS ay = 33 Saf m)! (l— 
@b od 3D Dont ye mi dem) 


( (3 II “TI IV 
(ft on (rab | 11) Soa On) (rac|12) OT (raa| 2) ret dr, + 


0 


P (28) 
“TI IV ist ve 
a fm (rac| 11) 01 (taal 171) 7 dry | er) on (Tab | 2) ata 
0 


0 
Finally we will remark that another important molecular integral, the three- 
center integral 


J Pa (1) Yr (1) Fe (1) dry (29) 
has been treated by means of the same method by Lowonrn [19], p. 41—43. 
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IV. Comments 


In this paper we have deduced a general potential formula, which has been 
applied to some two-centre, three-centre, and four-centre integrals of S-func- 
tions. These formulae, (23)—(28), are particularly suitable for computations by i 
means of numerical integrations if the start functions g(r) are given numerically, — 
either in the form of Hartree-Fock functions or in the form of tabulated ana- 
lytical functions. In general, the convergency of the infinite series in h is 
rather good; in practical applications to molecules and crystals, we have found — | 
it sufficient to take only the terms corresponding to h = 0, 1, 2, and.3 into — 
account.t 

Finally, we wish to point out that, in treating the cohesive energy of crystals 
or molecules having a high degree of symmetry, considerable simplifications can 
sometimes be performed by summing the contributions from the three-centre 
and four-centre integrals over suitable configurations of the nuclei involved. 
A very simple method of treating such symmetrical problems has been developed — Z 
by one of us (P.O. L.) and will be published elsewhere. 


Institute of Mechanics and Mathematical Physics, University of Uppsala, — 
Uppsala. 
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